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Abstract 

We present a two-dimensional lattice gas with anisotropic interactions which model 
the known properties of the surface reconstructions of CdTe and ZnSe. In contrast 
to an earlier publication JT2| the formation of anion dimers is considered. This alters 
the behaviour of the model considerably. We determine the phase diagram of this 
model by means of transfer matrix calculations and Monte Carlo simulations. We find 
qualitative agreement with the results of various experimental investigations. 

Keywords: Equilibrium thermodynamics and statistical mechanics, Monte Carlo 
simulations. Surface relaxation and reconstruction, Surface thermodynamics (includ- 
ing phase transitions), Cadmium telluride. Zinc selenide. Low index single crystal 
surfaces. 

1 Introduction 

Within the last years, potential technological applications of electronic devices based on II- 
VI semiconductors |jl| have inspired basic research concerning surfaces of these materials. 
In this context, various studies have addressed the properties of surface reconstructions. 
Experimental studies have investigated which reconstructions are present (2|, and how 
the reconstruction of the surface is influenced by parameters like temperature and particle 
flux in an MBE environment @, P . The majority of this work has been devoted to CdTe 
and ZnSe, where a fairly complete qualitative overview over the phase diagram has been 
gained. An overview over the properties of CdTe can be found in 0. On the other hand, 
there have been theoretical investigations of the reconstructions of CdTe [S] and ZnSe 



|10| using density functional theory. In these studies, knowledge about the chemical 
bonding of surface atoms and ground state energies of various reconstructions has been 
gained. 

Being based on quantum mechanics, density functional theory is believed to be exact 
apart from approximations made in the calculation. The computational burden of this 
method is comparatively high, which restricts its practical applicability to systems con- 
sisting of only a few atoms. Due to the periodicity of crystal surfaces this is not a severe 
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Figure 1: Panels (a),(b),(c): Sketches of the reconstructions of CdTe which are discussed 
in this work. The grey rectangles show the surface unit cells. The [110] axis is aligned 
horizontally, (a), (b): Cd-terminated reconstructions. Panel (a) shows the c(2 x 2)cd 
reconstruction, panel (b) the (2 x l)cd reconstruction. Panel (c) shows the (2 x l)Te 
reconstrution of surfaces terminated with a complete monolayer of Te. Panel (d) shows 
the attractive couplings in our model. 






restriction if one is interested in ground state properties of the system which however, are 
strictly relevant only at zero temperature. At higher temperature, the properties of the 
surface will be influenced by thermodynamic effects. These are particularly important if 
phase transitions between different reconstructions occur. Their theoretical investigation 
requires the study of systems with a large number of atoms which is beyond the scope 
of first principles methods or molecular dynamics simulations using realistic empirical 
potentials. 

Consequently, simplifying models are needed which preserve essential features of atomic 
interactions and can be investigated with moderate numerical effort. In many cases, 
two-dimensional lattice gases have been used successfully to model atoms adsorbed on a 
singular crystal surface or the terminating layer of such a crystal |11, ll^, |l^, 14]. In spite 



of the conceptual simplicity of such models the interplay of attractive and repulsive short 
range interactions can result in highly nontrivial critical behaviour and complex phase 
diagrams. In this paper, we will follow this approach to model the reconstructions of 
(001) surfaces of CdTe and ZnSe, our main focus being on CdTe. 

The outline of this paper is as follows: In section we will give a short review of the 
known facts about the reconstructions of the (001) surfaces of CdTe and ZnSe. In section 
^, we introduce a lattice gas model which considers the occupation of Cd sites and the 
dimerization of Te atoms and discuss its phase diagram. We conclude with a comparison 
of the features of our model with experimental results in section 0. 



2 Surface reconstructions of CdTe and ZnSe 

Both CdTe and ZnSe crystallize in the zinc-blende lattice. This lattice structure is com- 
posed of alternating layers of cations and anions which are parallel to the (001) surface, 
such that an ideal (001) surface would be terminated by a complete layer of one particle 
species. The positions of the atoms in one layer lie on a regular square lattice with its 
axes oriented in the [110] direction and the [110] direction. 

Under vacuum, the (001) surface of CdTe is Cd terminated. The surface is character- 
ized by vacancy structures where less than one half of the potential Cd sites in the top 
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layer are occupied |6[. This can be understood from simple quantum-mechanical con- 
siderations like the electron counting rule [15, 16], which states that a surface terminated 
by a complete layer of Cd is energetically unfavourable. At low temperature, one finds a 
c(2 X 2)cd reconstruction ||^, where Cd atoms and vacancies arrange in a checkerboard 
pattern (figure ||a). Frequently, a contribution of a (2 x l)cd arrangement can be found. 
In the (2 x l)cd structure, the Cd atoms arrange in rows along the [110] direction which 
alternate with rows of vacancies (figure |l|b). Density functional calculations |8| have shown 
that the surface energies of the two competing reconstructions are nearly degenerate and 
differ only by a small amount AE ~ 0.016ey per (1 x 1) surface unit cell. At a tempera- 
ture T ~ 570K, there is a phase transition [|2| above which the (2 x l)cd arrangement of 
the Cd atoms dominates the surface. An analysis of high resolution low energy electron 
diffraction (HRLEED) peaks [^] has shown, that there is a high degree of disorder in this 
high temperature phase. One finds elongated domains with a large correlation length of 
~ 375A in the [110] direction and a domain width as small as 22A in the [110] direction. 

If CdTe is exposed to an external Cd flux in an MBE chamber, the c(2 x 2)c(j recon- 
struction is stabilized at temperatures above the transition in vacuum. Under a Te flux, 
the surface is Te terminated with a (2 x l)Te reconstruction. At small Te fluxes, the surface 
is terminated by a complete monolayer of Te. The Te atoms on the surface form dimers 
which arrange in rows (figure ||c). At high Te flux and low temperature, one additional 
Te atom is incorporated into each dimer, such that the Te coverage of the surface is 1.5. 
The symmetry of this reconstruction is still (2 x 1), since the Te trimers tend to arrange 
in rows. A schematical phase diagram of the surface can be found in [^. 

Qualitatively, the properties of the reconstructions of ZnSe [17, ^] are quite similar 
to those of CdTe, where the Zn atoms are the counterparts of the Cd atoms and the Se 
atoms those of Te, respectively. There is one important exception: to date, no (2 x l)zn 
reconstructed surface has been found at high temperature. Density functional calculations 
[^, |l^ yield a higher energy difference AE f« O.OSeV per (1 x 1) surface unit cell between 
ideal c(2 x 2) and (2 x 1) reconstructed surfaces, which is approximately twice the value 
calculated for CdTe. As we will show, this greater energy difference might explain the 
different behaviour of CdTe and ZnSe. 



3 The lattice gas model 

The basic structure of our model of (001) surfaces is the same for different II- VI semicon- 
ductors. The differences between the materials are represented by the numerical values 
of parameters. In the following, for simplicity we will loosely speak of Cd and Te atoms 
instead of "cations" and "anions" without restricting ourselves to a modelling of CdTe 
only. 

In the lattice gas picture, a crystal lattice which is fixed in space is considered. Each 
lattice site is either occupied by an atom or empty. There are effective interactions be- 
tween atoms which include effects of surface strain and, at elevated temperature, lattice 
vibrations. Therefore, there is no simple mapping between ground state energies of the 
lattice gas model and surface energies determined from density functional theory. 

We model a flat (001) surface of CdTe i.e. we neglect the influence of steps on the 
reconstruction. This is a reasonable approximation if the typical distance between steps is 
much greater than the size of the unit cell of the reconstruction. In thermal equilibrium, 
this is fulfilled at temperatures well below the roughening transition. In this case, in the 
absence of bulk vacancies or other defects the crystal is uniquely described by the state 
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of the topmost monolayer which consists of one Cd half-layer and a Te half-layer. In our 

model, we consider such a monolayer with the Cd atoms on top. For simplicity, we assume 
the Tc layer to be fully occupied. This is not a severe restriction, since the removal of Te 
atoms will uncover Cd atoms in the layer below. Within the limit of a model of a flat 
surface, we simply do not distinguish whether a Cd terminated surface has been created 
by removing a Te layer or by depositing additional Cd atoms onto an intact Te layer. 

We use cartesian coordinates where the x-axis points in the [110] direction, the y-axis 
points in the [110] direction and the z-axis points in the [001] direction. The origin of 
the coordinate system is at the center of a Cd atom. Since we consider a flat surface, the 
z-coordinate of all atoms of one species is identical and will be omitted in the following. 
Measuring the lattice constant in appropriate units, the Cd atoms are at positions (x, y) 
with integer x,y. The Te atoms are at positions {x,y + 1/2), where x,y G Z. Tellurium 
dimers are created by the formation of a chemical bond between neighbouring Te atoms 
in y-direction ([110]). Since this is also the direction of the bonds of a surface Cd atom, 
a Tc dimer is formed by a pair of atoms, which might also be the binding partners of 
a Cd atom above them. This suggests a simple lattice gas representation both of the 
occupation of Cd sites and Te dimerization: We consider a rectangular array x of integers 
{•^'j}fj=:i- '^ij' ~ ^ represents a Cd atom at site which is bound to the Te atoms 

at sites — 1/2) and + 1/2). Xjj = 2 corresponds to a dimerization of Te atoms 

— 1/2) and + 1/2). Otherwise, Xjj = 0. 

In principle, one might also consider the formation of Te trimers by introducing a 
fourth state Xj j = 3. This corresponds to a trimer which consists of the Te atoms at 

— 1/2) and (i, j + 1/2) and an additional Te atom at (i,j). However, the formation 
of Te trimers plays a role only if the surface is exposed to a strong flux of pure Te at 
comparatively low temperature. Therefore, we have neglected this effect to reduce the 
number of parameters of our model and the numerical effort of its investigation. 

3.1 Interactions of atoms and dimers 

The detailed representation of the surface energies of a II- VI compound certainly would 
require long-range interactions and terms which depend on the simultaneous occupation 
of three or more sites. However, it is plausible to assume that the dominant contribution 
to the surface energy stems from pairwise interactions of atoms at short distances. In this 
section, we introduce a Hamiltonian which considers pairwise interactions between Xij on 
nearest neighbour sites and diagonal neighbour sites. These reflect the known properties 
of the reconstructions of CdTe and ZnSe. 

Due to the electron counting rule, surfaces with Cd coverages greater than 1/2 are un- 
stable while both a c(2 x 2)cd reconstruction and a (2 x l)cd reconstruction are permitted. 
This feature can be captured by introducing a hardcore repulsion between Cd atoms on 
neighbouring sites in the y-direction. In the x-direction, an attractive interaction favours 
the occupation of nearest neighbour pairs the strength of which is denoted by ex < 0. An 
attractive interaction ed <^ between diagonal neighbours stabilizes the c(2 x 2)cd recon- 
struction. These parameters are chosen such that the energy difference /S.E = |2erf — ea;|/2 
between these two reconstructions is small compared to the total surface energy per site. 

The electron counting rule favours Te dimerization, but forbids the formation of addi- 
tional bonds of dimerized Te atoms. Therefore, we forbid the simultaneous occupation of 
neighbour sites in the y-direction with dimers (formation of chains of Te-Te-bonds) and 
Cd atoms next to a dimer. These rules permit both a (2 x l)Te reconstruction, where the 
dimers arrange in rows and a c(2 x 2)Te reconstruction, where they arrange in a checker- 
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board pattern. Density functional calculations have shown that the surface energy 
of (2 X l)xo is significantly lower than that of a c(2 x 2)to reconstruction. This may be 
understood from the more efficient relaxation of surface strain in (2 x l)Te : Since the 
lattice is deformed in the same direction by dimers on neighbouring sites in x-direction, 
such an arrangement will be energetically favourable. We consider this fact by an attrac- 
tive nearest neigbour interaction ej between dimers in x-direction. Additionally, a dimer 
contributes a binding energy to the surface energy. Apart from the hardcore repulsion, 
we consider no interaction between Cd atoms and dimers. An overview over the couplings 
in our model can be found in figure |l]d. 
This yields the Hamiltonian 

L,N 

H = ^ exCijCi+ij + €dCij (cj+ij_i -I- Cj+ij+i) -I- €tdijdi+ij + Sbdij — fJ-Cij, (1) 

where we have introduced Cij = Sx^^^i and dij = j,2- The boundary conditions are 
assumed to be periodic. Cij represents the occupation of lattice sites with Cd atoms: 
Cij = 1 at the positions of Cd atoms and zero otherwise. Similarely, the positions of 
dimers are given by dij = 1. ^ is the effective Cadmium chemical potential, which 
includes the binding energy of surface Cd atoms to the substrate. The groundstate of the 
system at T = is determined by the chemical potential fi. For fi > fiQ = 2ed — et — eb, 
the surface configuration with minimal energy is c(2 x 2)cd • For ji < fiQ the groundstate 
is (2 X l)Te . 

In the following, we measure energy in dimensionless units which have been adjusted 
such that ed = —1- Additionally, we set A;^ = 1 such that temperature is measured in 
units of \ed\- 

3.2 Characterization of the surface configuration 

To characterize the surface reconstruction quantitatively, we evaluate the mean Cd cover- 
age pcd = the correlations 

measures the probability to find two diagonal neighbour sites which are simultaneously 
occupied by Cd atoms. This is a measure for the fraction of the surface which is covered 
by regions with a local c(2 x 2)c(j order. Similarely, C^^ measures the contribution of 
locally (2 x l)cd reconstructed regions in the system. 

The long range order of the Cd atoms is measured by the order parameters 

L,N 

^cf^ = I]^E^M^o^(^(^+J')) (4) 

^Cd"^^ is the staggered magnetization of a system of Ising variables {sij}fj^i, where 
Si J = 1 if a J = 1 and —1 otherwise. Large absolute values indicate a long range order 
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of the c(2 X 2)cd reconstruction. Its counterpart M^^^^^ measures the long range order of 
(2 X l)cd . 

Further, we introduce the mean Tellurium dimerization pu := {dij)-j and the dimer 
correlation Cf) := (dijdi^ij)- j, which characterizes the number of dimers which are in- 
corporated in locally (2 x 1) ordered areas. Their long range order is characterized by 

L,N 

3.3 Methods of investigation 

We have investigated this model by means of the transfer matrix method and Monte Carlo 



simulations. An introduction to the transfer matrix technique can be found e.g. in 21]. 
This method allows for a numerical calculation of the free energy of lattice systems with 
short-range interactions. In general, the system size must be finite in all directions but 
one. In our investigations, we choose the y-direction as the infinite direction. Since the 
computational effort increases exponentially fast with the system size L in the x-direction, 
only comparatively small L < 10 are feasible. However, this is sufficient to investigate the 
phase transitions of the system in the limit L ^ oo via finite size scaling techniques. To 



this end, we have applied the method suggested in |20| to determine the loci of first order 



phase transitions and to compute coverage discontinuities. Continuous transitions have 



been investigated by means of phenomenological renormalization group theory ||l^, 21 1 



Additionally, we have performed Monte Carlo simulations of our model. To obtain a 



reasonably fast equilibration, we have applied continuous time algorithms |22|. A more 
detailled description of these methods is presented in the appendix. We have simulated 
both the grand-canonical ensemble where pcd is controlled by p and the canonical ensemble 
where the number of Cd atoms in the system is fixed. In general, the results of the transfer 
matrix calculations and the Monte Carlo simulations are in good agreement (figure ^). 

3.4 Results 

To get insight into the typical behaviour of the model, we will first present a detailed 
investigation using the parameter set ed = €b = —1, ex = et = —1.9. The choice ex = —1.9 
reflects the fact that the energy difference between c(2 x 2)c(j and (2 x l)cd is small. 
The choice eb = —1 yields low concentrations of Te atoms which are neither dimerized 
nor bound to Cd atoms. This is consistent with experimental results. For simplicity, we 
have chosen typical example of the case where both couplings are of the same 

order of magnitude. However, since there is no next-nearest neighbour interaction between 
dimers, there is a significant difference in the surface energies of (2 x l)xe and c(2 x 2)Te, 
which is consistent with the results of [Hi- 



In section 3.4.4 we discuss the influence of the numerical values of the parameters on 
the phase diagram. We will show, that a smaller value of et affects the phase diagram 
qualitatively. In particular, this concerns the high temperature phases. 

3.4.1 Grand-canonical ensemble 

With the above parameter set, the Cd and the Te terminated groundstates are separated 
by fio = 0.9. In figure § the temperature dependent behaviour of the system at constant 
chemical potential is shown for p = 1 (figure ^,b) and fi = 0.8 (figure §c,d), which are 
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Figure 2: Simulations at constant chemical potential. Panels (a), (b) show the behaviour 

of the system with the couplings = e;, = — 1, e^; = = — 1.9 and a cadmium chemical 
potential /U = 1. With these parameters, the ground state of the system is c(2x2)cd • Panel 
(a) shows the Cd coverage pcd (+) and the correlations Cg^ (x), C^^ {*) and C'2 (Q). 

('2 X ('2x1) 

Panel (b) shows the mean absolute of the order parameters (©), M^^ (■) and 

Mjj ' (•). The symbols show data from a simulation run at a system size L = iV = 64. 

4 • 10^ • LN events have been performed both for equilibration and for measurement. The 
lines have been obtained by means of a transfer matrix calculation with a strip width 
L = 10. The data shown in panels (c), (d) have been obtained at = 0.8, where the 
ground state is (2 x l)Te • AH other parameters are identical to those used in (a), (b). 
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examples for both cases. At /x = 1, the ordered c(2 x 2)cd phase at low T manifests itself 
in a high Cd coverage pcd and values of the correlation C^^ and the order parameter 

^Cd"^^ close to 0.5. In figure ^ a surface snapshot in this phase at T = 0.71 is shown. 
As expected, the Cd atoms arrange preferentially in a checkerboard configuration. At 
a temperature T = 0.83, there is a first order phase transition to a disordered phase. 

(2x2) rl 

This is indicated by a decrease of the order parameter and the correlation C^^. 

Simultaneously, C^^ starts to increase. At T = 1, both lines cross such that the high 
temperature behaviour of the system is dominated by a local (2 x 1) ordering of the Cd 
atoms. The simulation data shown in figure 0a, b have been taken from two independent 
simulation runs. The system was initialized with a perfect c(2 x '2)q^ configuration at the 
lowest temperature investigated. Then, as successively higher temperatures were imposed, 
the surface configuration was kept as initial state of the next simulation. This reflects in 

(2 X 2) 

a small hysteresis effect in due to the first order nature of the phase transition. 

At = 0.8, we obtain a completely different behaviour. Since the groundstate is a 
(2 X l)xe reconstruction, at small T we measure low Cd coverages and values of C|j and 
close to 0.5. The most frequent thermal excitations are Cd adatoms, the density 
pcd of which increases with T. Figure ^ shows a surface snapshot at T = 0.93. The 
Cd atoms preferentially arrange in rows such that ^ ^cd- These rows adapt to the 
structure of the (2 x l)Te reconstruction. This yields nonzero values of the order parameter 
M^^^\ which indicate a global ordering of Cd atoms in a (2 x 1) arrangement. However, 
the interactions between the Cd atoms themselves are insufficient to stabilize this global 
order. Instead, it is purely induced by the environment of the Te dimers. 

At T = 1.67, there is a first order phase transition above which the system is in a 
disordered phase similar to that found at ^ = 1. Remarkably, we observe high pcd ~ 0.2 
at temperatures slightly below the phase transition. At the phase transition, pcd jumps 
to an even higher value. Above the transition it decreases slightly with T. 

3.4.2 Phase diagram 

Figure ^,b shows the phase diagram of the model, which has been extrapolated from 
transfer matrix calculations with strip widths L of 6, 8 and 10 lattice constants. In figure 
^b, the lines of phase transitions in the p-T plane have been plotted. At low temperature, 
the system is either in an ordered (2 x l)Te (1) or an ordered c(2 x 2)cd phase (2). The line 
of phase transition between these phases starts at at zero temperature and p = po, where 
the energies of both reconstructions are degenerate. For < T < Tt = 0.84, it remains at 
the same chemical potential po, apart from small numerical uncertainties of extrapolation. 
In consequence, a phase transition between a c(2 x 2)cd and a (2 x l)Te reconstruction 
cannot be observed at constant chemical potential. The disordered phase (3) exists for 
T > Tj. At the point (T = Tt,p = po) five phases coexist: the disordered phase, and 
c(2 X 2)cd and (2 x l)xe reconstructed phases in two sublattices, corresponding to positive 
and negative values of the order parameters M^J''^^ and AI^^^^ , respectively. The (2 x l)Te 
reconstructed phase (1) exists only at temperatures below a critical temperature = et- 
At T^, the line of the phase transition to the disordered phase (3) diverges to = — oo. On 
the contrary, the c(2 x 2)cd reconstructed phase (2) may exist at arbitrary temperature, if 
the Cd chemical potential is large enough. The dashed line in figure ^b shows the chemical 
potential, at which the correlations Cq^ and Cq^ in the disordered phase are equal. For 
smaller p, the local ordering of Cd atoms is dominated by a (2 x 1) arrangement, while 
for larger p they prefer a local c(2 x 2) ordering. 
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Figure 3: Phase diagram with parameters = e;, = —1, e^^ = et = —1.9. The sohd hnes 
in panel (a) show the hnes of phase transitions as functions of pcd and T, in panel (b) 
they are shown in the [i-T plane. Note the offset in the temperature axes. In region 
(1), the system is in a homogeneously ordered phase with a (2 x l)Te reconstruction, in 
region (2) it is homogeneously ordered and c(2 x 2)cd reconstructed. Region (3) is the 
disordered phase. On the left side of the dashed line, the Cd atoms show preferentially a 
local (2x1) ordering (Cg,^ > C%^, while on its right side a c(2 x 2) ordering dominates. 
Due to the coverage discontinuity between (1) and (2), (3), there is a coexistence regime 
where regions with high and low pcd coexist (4). The symbols show lines of constant 
[I (Same data as in Figure |2|a,c) Panels (c)-(f) show typical surface snapshots, (c)-(e) 
correspond to grand-canonical simulations. The snapshots show surface configurations 
after 8 ■ 10^ • LN events, (c): T = 0.93, /u = 0.8 (Phase 1), (d): T = 0.71, p = 1 (Phase 2), 
(e)T = 1.14, = 1 (Phase 3). (f) displays a surface configuration after 2 ■ 10^ • LN events 
in a canonical simulation at T = 0.73, pcd = 0.35 (Region 4, Coexistence of Phases 1 and 
2). All snapshots show sections of 32 x 32 lattice constants of systems of size L = N = 64. 
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Figure |3|b shows the phase diagram in the pcd-T plane. At the transition between 
the phases (2) and (3), pcd varies continuously. The phase transition from phase (2) to 
phase (1) occurs at a temperature independent chemical potential pQ. Therefore, there is 
no coverage discontinuity if temperature is varied at a constant chemical potential fi > fiQ 
where the groundstate is a c(2 x 2)cd reconstruction. 

On the contrary, there is a coverage discontinuity at the phase transition between 
the ordered phases (1) and (2) and at the transition between phases (1) and (3). These 
discontinuities yield a coexistence regime (4). Here, for T < Tt, c(2 x 2)cd and (2 x l)Te 
reconstructed phases coexist, while for T > Tt the (2 x l)Te reconstructed phase coexists 
with the disordered phase. Figure ^ shows a surface snapshot from a simulation which 
was performed at a constant Cd density pcd = 0.35 and a temperature T = 0.73. This is 
a typical surface configuration in the regime where the ordered phases (1) and (2) coexist. 
The system is separated in two phases, one with a high pcd and a c(2 x 2)cd reconstruction, 
and another one with a (2 x l)Te reconstruction and a low concentration of Cd adatoms. 
The local values of pcd in both regions are given by the left and the right boundary of 
the coexistence regime. At low temperature one obtains pcd ~ 0.5 in the c(2 x 2)cd phase 
and pcd ~ in the (2 x l)Te phase. At temperatures T ^ 0.6, pcd in the (2 x l)Te phase 
increases strongly with T and obtains its maximal value p\j2°^^ ~ ^-^^ at T = 0.98. At 
even higher temperature it decreases with T and becomes zero at T^. On the contrary, 
pcd in the c(2 x 2)cd phase remains high for T < Tt. At T > T^, the Cd-rich phase is 
disordered. Then, pcd at the right boundary of the coexistence regime decreases with T. 
At TI, it becomes zero and the coexistence regime disappears. The dashed line in figure |3|a 
marks the values of pcd at which C^^ = C^^. For smaller coverages, C^^ > C^^ such that 
the local ordering of the Cd atoms is dominated by a (2 x 1) arrangement, while for greater 
coverages, C^^ < C^^. The values of pcd which have been measured in the simulations 
shown in figure |2| are plotted in the phase diagram as examples of lines of constant chemical 
potential. The loci of the phase transitions and the point where = Cq^ (at p = 1) 
are in good agreement with the results of the transfer matrix extrapolation. 



3.4.3 Canonical ensemble 

The temperature dependent behaviour of the model at a constant Cd adatom density 
PCd = 0.35 is shown in figure This is an example which shows the typical behaviour of 
the model under the conditions of phase separation. At low temperature, most of the Cd 
atoms are concentrated in a Cd rich phase with a c(2 x 2)cd reconstruction. The remaining 
area of the system is covered with a (2 x l)xe reconstructed phase. Since both phases are 

(2x2) f2xl) 

long-range ordered, we obtain large values of the order parameters M^^ and Mj^ and 
the corresponding correlations. The fraction of Cd atoms which occupy sites in the Te rich 
phase yields nonzero values of the order parameter M^^ . As temperature increases, more 
and more Cd atoms pass into the Te rich phase. This yields an increase in Cq^ and 
and a decrease of C^^ and Mj^^"^^ due to the (2 x 1) ordering of the Cd atoms in the Te rich 
phase. At the temperature Tj, the Cd rich phase undergoes the order-disorder transition 

(2 X 2) 

where M^j^^ drops to zero. At this phase transition, the fraction of Cd atoms which are 
incorporated in locally (2 x 1) ordered configurations in the Cd rich phase increases. Thus, 
there are two independent effects which lead to a dominance of C^^ over at higher 
temperature: The specific shape of the left boundary of the coexistence regime and the 
order-disorder transition of the Cd rich phase. At a temperature T^, the system is leaving 
the coexistence region into phase (3). At this phase transition, the separation between a 
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Figure 4: Results of a simulation of a 64 x 64 system at conserved pcd = 0.35 and couplings 
= Cft = — 1, Ca; = = —1-9. 10^ • LN events have been performed both for equilibration 
and for data sampling. Panel (a) shows the correlations Cq^ (x), C^^ (*) and (□). 

Panel (b) shows the mean absolute of the order parameters (©), M^^ (■) and 

M^^^-* (•). The inset shows standard deviations of order parameters. The meaning of 
the symbols is the same as in the large picture. 



Te rich and a Cd rich phase disappears, such that at high temperature the system is in 
a homogeneous, disordered state. The order parameters M^^^^ and M^^^^ become zero 

(2x1) 

and the local correlation between dimers, Cjj ' , decreases. In the simulations we have 
determined the critical temperatures of the phase transitions from the standard deviations 
of the order parameters M^^"^^ and M^^^\ Mq^"^^ is peaked at Tt, where the long range 
order of the Cd atoms is lost, while is peaked at T^, where the system leaves the 

coexistence regime and the dimers lose their long range order. We obtain Tt = 0.79 it 0.2 
and = 1.0 it 0.2, where the uncertainty is due to the temperature spacing between the 
single simulations. These values are systematically lower than the theoretical results of the 
transfer matrix extrapolation {Tt = 0.84, T^ = 1.11). However, the theoretical results are 
valid in the limit of an infinite system size, where the free energy of the phase boundary 
can be neglected compared to that of the bulk of the phases. We have verified that the 
systematic deviation between theory and simulations decreases with the system size. 



3.4.4 The influence of the parameter set on the phase diagram 

Finally, we discuss the influence of the parameter set on the phase diagram. The main 
effect of a variation of the binding energy of the Te dimers e;, is to shift the chemical 
potential at which the phase transitions occur. The influence on the shape of the phase 
diagram is small as long as e;, is sufficiently large such that the density of Te atoms which 
are neither dimcrizcd nor bound to Cd atoms is small. Due to the electron counting rule, 
a state where a large fraction of the Te atoms remains unbound should be irrelevant for 
CdTe. 

A smaller energy difference A.E between a perfect c(2 x 2)cd and a perfect (2 x l)cd 
reconstruction increases the tendency of Cd atoms to arrange in a 2x 1 order. If = —1.95 
and all other parameters are identical to our standard parameter set, the dashed line 
where C^^ = C^^ is shifted to a higher pcd ~ 0.45. Additionally, we obtain a higher 
value of the maximal Cd coverage in the Te rich phase Pq^"'^ = 0.28, which is achieved 
at a lower temperature T = 0.9. The temperature Tt = 0.7 above which the disordered 
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Figure 5: Phase diagram with the parameter set ed = eb = —1, = —1.6, = —1.2. 
Compared to the phase diagram shown in figure ^, the energy difference between the 
c(2 X 2)cd and (2 x l)cd reconstruction is greater and the interaction between dimers is 
weaker. With these parameters, the transition between phases (1) and (3) is continuous. 
Panel (a) shows the phase diagram in the pcd-T plane, in panel (b) it is shown in the p-T 
plane. 

phase is stable is also slightly lower. In general, a greater value of le^^l increases both 
the tendency of the Cd atoms to arrange in a (2 x 1) order in the disordered phase and 
the concentration of Cd atoms in the (2 x l)xe reconstructed phase. The temperature 
of the order-disorder transition in the Cd rich phase is lowered. Conversely, a greater 
energy difference between perfect c(2 x 2)cd and (2 x l)cd reconstructed phases alters the 
temperature of this transition and leads to a preferential c(2 x 2) arrangement of the Cd 
atoms in the disordered phase. 

The main effect of a variation of the interaction et between Te dimers is to change the 
properties of the transition between the phases (1) and (3). As long as et is sufficently 
high, the phase diagram is qualitatively similar to that shown in figure ^. In this case, the 
main effect of a variation of et is to shift the temperature ~ et where the (2 x l)Te phase 
vanishes. However, at low et this transition becomes a continuous phase transition without 
any coverage discontinuity. As an example, in figure ^ a phase diagram with et = -1.2, 
Cfe = = — 1 and ex = —1.6 is shown. The coexistence regime (4) vanishes at Tt = 0.9 
such that there is no phase separation between the long-range ordered (2 x 1)tc phase and 
the disordered phase. This parameter set yields a comparatively high Ai?. Therefore, the 
line where C^^ = C^^ is at relatively low Cd coverages. In contrast to the situation at 
large et, the chemical potential at which the transition from phase (1) to phase (2) occurs 
is not independent of temperature. Instead, the transition is at slightly greater ^ at higher 
temperature. In consequence, there is a small range 0.2 < fi < 0.27 where there is a phase 
transition from phase (2) to phase (1) if T is increased at constant chemical potential. 

4 Comparison with experimental results 

The results presented in section |3.4| suggest an interpretation of the experimentally ob- 
served crossover from a c(2 x 2)cd reconstruction to a (2 x l)cd reconstruction as an 
accompanying effect of an order-disorder phase transition. At low temperature, there is 
a long-range ordered Cd-rich phase with a c(2 x 2)cd reconstruction. At a critical tem- 
perature, the Cd atoms lose their long-range order and arrange preferentially in a (2 x 1) 
pattern. This picture is consistent with the experimental observation of small domains in 
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the (2 X l)cd reconstruction Q which indicate a high degree of disorder. 

Strictly speaking, a CdTe surface under vacuum is not in thermal equilibrium. At 
the temperature of the c(2 x 2)cd - (2 x l)cd transition, sublimation plays an important 
role. However, in a previous publication we have shown that the basic features of a 
simplified version of our model which neglects Te dimerization [^] are preserved under the 
conditions of step flow sublimation. This is the dominant sublimation mechanism for CdTe 
(001) 1^^. Sublimation is slow enough to permit a local equilibration of the terminating 
layer which justifies the application of equilibrium thermodynamics. However, the Cd 
coverage pcd is determined by the sublimation process. Therefore, within the limit of an 
equilibrium model it is not possible to calculate the path in the phase diagram that CdTe 
follows. 

Auger measurements [^] yield pcd ~ 0.35 at the transition temperature in vacuum. 
This suggests that the system is in the coexistence regime for a wide range of temperatures. 
Then, the behaviour of CdTe in vacuum should be similar to our results at constant 
PCd = 0.35 (section 3.4.3 ). Electron diffraction techniques investigate large regions on 
the surface and hence yield averages over all coexisting phases. Consequently, it should 
be reasonable to compare results of these experiments with quantities which are averaged 
over the whole system in our canonical simulations. As discussed in section 3.4.3, there are 
two independent effects which lead to a clear dominance of the arrangement of Cd atoms 
in rows: the order-disorder transition in the Cd-rich phase and the fact that an increasing 
fraction of the Cd atoms is dissolved in the Te-rich phase. Electrons are diffracted both 
from the Cd atoms and the Te dimers which have (2 x 1) order. The superposition of these 
effects should yield the pronounced (2 x 1) diffraction peaks which have been observed in 
I, I- 

In the c(2 x 2)c(j reconstructed phase, at temperatures well below the phase transition 
one frequently finds collective thermal excitations of adatoms, where a row of several Cd 
atoms is shifted by one lattice constant in the y-direction. Due to the repulsion between Cd 
atoms in this direction, this is possible only if one Cd atom per excitation is missing. An 
example is shown in figure |3|d. This effect has been observed experimentally by Seehover 
et. al. 1 25] at room temperature using STM microscopy. There is a striking similarity 
between figure 3 in |25] and figure 

The phase diagram of our model explains the properties of the reconstructions of the 
CdTe (001) surface under an external particle flux. If the Cd coverage of the surface is 
increased by deposition of Cd, the state of the system moves into regions of the phase 
diagram where the Cd atoms arrange preferentially in a c(2 x 2) pattern. Depending on 
temperature, this is either the ordered c(2 x 2)cd phase (2) or the region of the disorded 
phase (3) on the right side of the dashed line where C^,^ > C^^. Indeed, experiments 
[Q, have shown that an external Cd flux restores the c(2 x 2)q^ reconstruction at high 
temperatures where a (2 x 1) order is found under vacuum conditions. Clearly, a strong Te 
flux induces the formation of a long-range ordered (2 x l)Te phase at temperatures below 

c 

On ZnSe, a Zn terminated (2 x 1) reconstruction has not been observed yet. Our 
model offers an explanation of this fact which is consistent with the results of density 
functional theory |l^, |8|. These calculations have shown, that the difference in the 
surface energies per (1 x 1) surface unit cell between perfect cation terminated c(2 x 2) 
and (2 x 1) reconstructions in ZnSe is significantly greater than in CdTe (0.03eV versus 
O.OOSeV). In our model, this greater energy difference corresponds to a smaller value of 
\ex\ which shifts the line where = to smaller coverages. Consequently, a (2 x 2)zn 
arrangement dominates in a much wider range of coverages. 
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Wolframm et. al. have measured the locus of the transition beween the (2 x 2)zn 
reconstruction and the (2 x l)se reconstruction as a function of temperature and the com- 
position of an external particle flux by means of reflection high energy electron diffraction 
(RHEED). They find that at higher temperature a greater Se flux is required to obtain 
a (2 X 1) diffraction pattern. At temperatures above 450 °C no (2 x l)se reconstruction 
could be observed even under extremely Se-rich conditions. This is reminiscent of our 
observation that the anion-rich phase (1) vanishes at a temperature T):. 

Unfortunately, the available experimental data are insufficient for a systematic fit of 
the model parameters. However, some rough estimates show that at least the orders of 
magnitude are reasonable. As discussed above, the c(2 x 2)cd -(2 x 1)qa transition of CdTe 
in vacuum should be at a temperature close to Tt. Identifying this with the experimental 
value of 570 K, we obtain the value of our energy unit \ed\ ~ 0.06eV in physical units. 
This yields a value /S.E ~ O.OOSeV for the difference in the surface energies of c(2 x 2)cd 
and (2 x l)cd i which is about 1/2 of the value of O.OOSeV which has been obtained by 
means of DFT calculations. This shows at least that the qualitative agreement between 
experiments and our model has been obtained in a physically reasonable region of the 
parameter space. 

In our model, phase (1) vanishes at a temperature ~ |ef|. Identifying this with 
the value of 450°C which has been measured in experiments on ZnSe, we obtain that 
et ~ O.OGeV. This is the same order of magnitude as our estimate of \ed\ in CdTe. 

These considerations suggest that it should be possible to obtain quantitative agree- 
ment between experiments and our model both for CdTe and ZnSe with values of the 
model parameters on the order of magnitude of a few ten meV. 
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Appendix: Algorithms for continuous time Monte Carlo sim- 
ulations 

We have simulated our model using continuous time Monte Carlo techniques which greatly 
improve the computational speed compared to simpler, Metropolis like algorithms ||2^ . A 
Markov chain of states s of the system is constructed with transition rates Vg^s' between 
states s and s'. This dynamics converges to a Gibbs distribution if every possible state of 
the system can be reached within a finite number of transition events (ergodicity) and the 
rates fulfil a detailed balance condition r^-^s'/rs'^s = exp(— (i7(s') — H{s))/T). In every 
time step, one event is performed which is selected randomly with probability Ps^s' = 
r s^s' / R{s), where R{s) = Yls'''^s-^s' is the sum of transition rates of all possible events 
in state s. The physical time interval At{s) = 1/R{s) which passes between subsequent 
events depends on the state of the system. In the calculation of thermal averages, the 
observables measured in state s have to be weighted with At{s). This algorithm requires 
the knowledge of the rates of all possible events in the current state of the system. If the 
configuration of the system is changed only locally in an event, both the selection of events 
and the updates of the rates of events can be done in O (log TV) CPU steps using a binary 
search tree. Here, J\f is the number of possible events. 



14 



The application to grand-canonical simulations is straightforward. There are two pos- 
sible events per lattice site which change the state of the site: Xij ^ (xij + 1) mod 3 and 
Xij <— {xij + 2) mod 3. These simulations can be done with Metropolis-type rates rg^s' = 
min{l, exp[— (i?(s') — H{s)/T]} or symmetrical rates rs^s' = exp[— (if(s') — H{s))/{2T)]. 
Both possibilities yield identical results. 

A canonical simulation requires an algorithm where the number of particles is fixed. 
In this case, in every event one particle jumps to a different site. We choose a nonlocal 
dynamics where the range of particle jumps is unlimited. This yields considerably faster 
equilibration compared to a Kawasaki dynamics with nearest neighbour diffusion only. 
For simplicity, we permit only jumps to a site where the binding energy of the particle is 
independent of the state of its initial site, i.e. we forbid jumps to nearest and next nearest 
neighbour sites. If the particle jumps from site i to site j, the energy difference between 
the final and the initial state is AH = AHj — AiJj, where AHx is the energy difference 
of the system with site x occupied and empty. The rates 

n^j = exp [{AHi - AHj) / (2T)] (7) 

fulfil the detailed balance condition. Then, the probability for a jump from site i to site j 
factorizes, i. e. 

Pi-^j = Pi ■ Pf ^'^^^^ Px = 

Here, we have introduced the rate for deposition of a particle at site x (r+) and for removal 
of a particle at site x (r^). = exp[—{AHx) / {2T)] if site x is empty and zero otherwise. 
Conversely, r~ = cxp[AHx/ {2T)] on occupied sites and zero on empty sites. Due to this 
factorization property we can proceed in two steps: In the first step, we select the site i 
from which the particle starts with probability . Then, we select the site j where the 
particle is landing with probability p^. If the distance between site i and site j is > 2 the 
particle is moved. Otherwise, the event is rejected and the system remains unchanged. 
Since the number of rejected events is small on large systems, the loss of speed can be 
neglected. 

However, in our model only the number of Cd atoms is fixed, while the number of Te 
dimers may change. This requires a more elaborate algorithm which uses both canonical 
and grand-canonical techniques. A Cd atom may jump to any site which is not occupied 
by a Cd atom. If the arrival site is occupied by a dimer, the dimer is destroyed. The 
removal of the Cd atom at the starting site creates a pair of Te atoms. We consider 
both the case where these Te atoms dimerize immediately and the case where they remain 
unbound. Additionally, Te dimers may break up and dimers may be created on empty 
sites. In the following, these processes will be denoted as "dimer flips" . In each timestep, 
a random number is drawn to decide whether a Cd jump with immediate formation of a 
dimer at the starting site, a Cd jump without dimerization or a dimer flip will occur. The 
probabilities are proportional to the sums of the rates of the processes in each group. A 
L X N system contains 0{LN) Cd atoms and Te dimers. Thus, there are 0{LN) dimer 
flip events. Since there are 0{LN) possible arrival sites for each Cd atom, there are 
0{L?'N'^) Cd jump events. To keep the ratio of dimer flips and Cd jumps independent of 
the system size, the dimer flips have been weighted with a prefactor LN . Then, one event 
in the selected group is perfomed. For dimer flips, this is done with the grand-canonical 
algorithm while Cd jumps are performed by the canonical two-step algorithm. 
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